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The theory of controUed quantum open systems describes quantum systems interacting with 
quantum environments and influenced by external forces varying according to given algorithms. It 
is aimed, for instance, to model quantum devices which can find applications in the future technology 
based on quantum information processing. One of the main problems making difficult the practical 
implementations of quantum information theory is the fragility of quantum states under external 
perturbations. The aim of this note is to present the relevant results concerning ergodic properties of 
open quantum systems which are useful for the optimization of quantum devices and noise (errors) 
reduction. In particular we present mathematical characterization of the so-called " decoherence-free 
subspaces" for discrete and continuous-time quantum dynamical semigroups in terms of C* -algebras 
and group representations. We analyze the non-Markovian models also, presenting the formulas for 
errors in the Born approximation. The obtained results are used to discuss the proposed different 
strategies of error reduction. 



There exists a widespread believe that sooner or later the technology of the future will employ typical quantum 
features of microscopic or mesoscopic systems like entanglement or superpositions of quantum states. The most 
prominent examples are the rapid development of the quantum information and computation theory [1,2] and the 
corresponding progress in experimental techniques. The later includes realization of quantum teleportation [3], mon- 
itoring of dissipation and decoherence for mesoscopic quantum systems [4], construction of quantum gates [5], etc. 

The theory and experiment show that entanglement of separated systems and superpositions of well distinguishable 
states are very sensitive to perturbations by environment. The interaction with environment introduces noise (errors) 
which should be taken into account by designing any quantum devices and their operations. The appropriate theory 
has been developed in the 70-ties and 80-ties as the Quantum Theory of Open Systems (see [6-11] for the survey 
of this theory). The basic idea of this theory is the decomposition of the Universe into the open system S and the 
environment (reservoir , heat bath) R with their corresponding Hilbert spaces T-Ls and Tiji respectively. Assuming 
that the interaction between S and R is weak we can always write, within a reasonable approximation, that for the 
initial moment to = the initial state of the total system is given by the product formula 



where p is an arbitrarily prepared initial state of S and ljr is a stable reference state of R. Our interest is concentrated 
on the time evolution of S for times t>Q given by the reduced density matrix (we omit subscript S) 



where {Ut\t € R} is a family of unitary operators acting on JisR ~ Ti-s ® "Hr and describing the reversible (Hamilto- 
nian) evolution of the total system. Tr^j is a partial trace over Ti.R. The dynamical maps {Af ; t > 0} , their properties 
and different approximation schemes for them (e.g. Markovian approximations) are the main topics of the Quantum 
Theory of Open Systems. 

In order to deal with the problems posed by quantum information processing and quantum computing one has 
to extend the standard approach introducing, beside the open system S and the reservoir R, the third element - 
the controler C. The role of C is to prepare an initial state , to measure a final one and to control, as far as 
it is possible, the time evolution of S by means of the time-dependent Hamiltonian Hs{t). The presence of the 
time-dependent Hamiltonian makes, for example, the application of Markovian approximation often impossible and 
demands new mathematical techniques. Moreover, the physical structure of C, the resources used for state preparation 
and measurements must be taken into account in the estimation of efficiency of quantum devices. We shall call such 
systems Controlled Quantum Open Systems (CQOP). 

One of the main problems in the theory of CQOP's is to reduce errors due to the interaction of a quantum device 
with an environment. At least three elements of the presented scheme can be optimized: the choice of the subspace in 
Hs which supports the initial states of S, the trajectory {Hs{s); < s <t} which leads to a given unitary operation 
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on Tis and the structure of the couphng of S to R. In the recent hterature one can find a number of proposals of such 
optimization schemes [12-14]. However, some of these ideas are essentiahy based on the independent derivations of 
the ear her results, some other employ oversimplified phenomenological models which should be scrutinized from the 
point of view of a more fundamental theory. Therefore, it seems to be useful to collect, improve and generalize the 
relevant results concerning mainly the ergodic properties of quantum open systems and discuss their consequences for 
the optimization schemes mentioned above. 

In the Section 2 wc present mathematical characterization of the subalgebras of observables which evolve according 
to reversible dynamics for a given irreversible quantum dynamical map or quantum dynamical semigroup. Such 
subalgebras correspond to " decoherence-free subspaces" in the language of modern quantum information theory. 
The general ideas arc illustrated by examples of collective dissipation models in SectionS. Section 4 is devoted to 
the derivation and discussion of the formula describing errors in CQOP in the non-Markovian weak coupling Born 
approximation. 

II. ERGODIC PROPERTIES OF QUANTUM DYNAMICAL MAPS AND SEMIGROUPS 

In this Section we restrict ourselves in most cases to quantum open systems with finite-dimensional Hilbert space 
Hs = C". We discuss the irreversible dynamics either in terms of the single dynamical map F and its iterations 
r'^, fc = 0, 1, 2, ... or the dynamical semigroup {Tt, t > 0} with the composition law 

TtTs=Tt+s ,t,s>0 , To = l . (3) 

Both r and Tt describe the dynamics in Heisenberg picture i.e. they act on a matrix algebra (C*-algebra) M.n 
representing quantum observables of the system. In the following all results independent on the semigroup property 
(3) are formulated in terms of a single dynamical map T. The Schrodinger picture evolution is given by the dual map 
r* acting on M.n which is treated now as the smallest complex linear space containing density matrices (quantum 
states) of the system and satisfying 

Tr(Ar(B)) = Tr(r*(A)S) . (4) 
Dynamical maps must preserve positivity and normalization i.e. 

T{A*A) > , T*{B*B) > , 

r(l) = 1 , TtT*{B) = Tv{B) (5) 

for ah A,B € M„. 

A. Complete positivity 

Positivity preserving condition (5) is generally to weak for a mathematically and physically consistent theory. The 
minimal condition which must be imposed on dynamical maps to allow meaningful construction of joint dynamics for 
noninteracting composed systems is complete positivity. It means that for any d= 1,2, 3, F (g) 7^ is positive,where 
Id is an identity map acting on d x d matrices (i.e. trivial dynamical map on d-lcvel quantum system). Any reduced 
dynamical map obtained from the formula (2) is completely positive and any completely positive dynamical map can 
be obtained from the reduced dynamics scheme. On the mathematical side completely positive maps on operator 
algebras were studied already in the 50-tics and the celebrated Stincspring representation [15] leads to a general form 
of completely positive dynamical map called often Kraus decomposition [16] 

t{a)=y;^w:aw^ (6) 

a 

where Wa are bounded operators on H satisfying J2a^a^a = 1- The decomposition (6) is highly nonunique, in 
particular the sum over a can be replaced by an integral. If H is n-dimensional then complete positivity is equivalent 
to n-positivity and one can always find Kraus decomposition in terms of at most terms. 
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From now on by the quantum dynamical semigroup (in the Heisenberg picture) we mean a a one-parameter strongly 
continuous semigroup Tt,t>0 of completely positive unity preserving maps. For finite dimensional case we have 
always 

At = Tt{A) = e'\A) , j^At = L{At) (7) 
with the generator of the standard Lindblad, Gorini, Kossakowski and Sudarshan form [17] 

L{A) = t[H,A]+j2v;AVj - lij^y/yj'^} 

= ^[H,A] + ^j2i^;[A,VJ] + [V;,A]VJ) . (8) 

The choice of the Hamiltonian H = H* and the operators (matrices) Vj is again not unique and the sum over j can be 
replaced by an int(^gral. In our finite dimensional case one can always find at most p = — 1 siich Vj-a. To simplify 
the notation we put always h=l and fcs = 1 to have the same units for energy, frequency and temperature. 

The irreversible character of a dynamical map or semigroup is expressed in terms of the following inequalities valid 
for any A e Mn 

T{A*A) > r(^*)r(^) - - Kadison inequality (9) 

L{A*A) > L{A*)A + A*L{A) - - dissipativity. (10) 

Indeed, puting equality in (9) we obtain in our case the reversible dynamical map A i-^ U*AU with a unitary U 
while the equality in (10) implies the Hamiltonian derivation structure L{A) = i[H,A]. Both inequalities are the 
consequences of 2-positivity condition obviously satisfied for any completely positive dynamical map and proved very 
useful in the analysis of ergodic properties. 

B. C*-algebras and group representations 

In this Section we briefly review the basic mathematical facts about finite dimensional C*-algebras and represen- 
tation theory of compact groups which will be used later on [18]. 

The full matrix algebra is an example of C*- algebra i.e. a complex linear space with adjoint operation (hermitian 
conjugation) A i— » A*, multiplication AB and complete with respect to the norm satisfying = ||^||^. We are 

interested in finite-dimensional C*-algebras which can always be seen as C*-subalgebras of full matrix algebras i.e. 
linear subspaces of M.n, closed with respect to multiplication of matrices and hermitian conjugation and containing 
the unit matrix 1. 

Any such C*-subalgebra Af can be represented as the algebra of linear maps on the finite dimensional Hilbert space 
in the following way. The Hilbert space is decomposed as 

m m 

C" = C"^- O , Ujdj =n (11) 

where nj is the dimension of the irreducible representation of J\f and dj its multiplicity. Then any element A G J\f 
can be written in the form 

m 

A = ^Aj<^U, (12) 

where Aj e and 1^^. is a unit matrix acting on C*^ . 

In a general algebraic framework the physical symmetry or reversible dynamical map described in the Heisenberg 
picture is given in terms of an automorphism of the corresponding C*-algebra, i.e. a linear map A i— > U{A) satisfying 
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U{A*) ~ U{A)* , U{AB) = U{A)U{B) , U{1) = 1 and reversible. For an automorphism of a finite dimensional algebra 
M C Mn reversibility follows from the former conditions and moreover there exists a unitary matrix U G Mn such 
that U{A) = U*AU. The group of all automorphisms of Af is denoted by Aut{M). A subgroup AutQ{M) C Aut{M) 
contains all inner automorphisms given in terms of imitarics U G M. For our finite dimensional case AuIqIM) 
coincides with the path connected component of identity of Aut{N). In order to illustrate the fact, that even in the 
finite dimensional case not all automorphisms are inner, consider the abelian subalgebra of M.2 of diagonal matrices 

and the "flip" map ^b) ^ i^Q a) ' ^^^^ automorphism can be writen as A ^ cFxAgx, where the (unitary) 

Pauli matrix is neither diagonal nor cannot be continuosly transformed into identity. 

Wc denote by A' = {X G Mn\ [X,A\ — 0} the commutant of the subset A C M.n- If A is sclf-adjoint, i.e. 
A ^ A ^ A* ^ A, then the commutant A' is a unital C*-subalgebra of For a subalgebra N we have (A/"')' = N ■ 
Obviously \i A ^ B then B' C A' . Denoting by Alg{A) the smallest C*- algebra containing A we have A' = Alg{A)' . 
For an algebra U = ©j^i Mn, ® Idj (see (12)) its commutant M' = ©^=i In, ® Mid, ■ 

The theory of group representations provides natural examples of finite dimensional C*-algcbras. Take a finite- 
dimensional unitary representation TZ : G ^-^ Mn', 9 ^ R{9) of the compact group G (discrete or Lie group). Such 
representation can be always decomposed into irreducible representations TZj according to (11). Any subspace C"' is 
a carrier of an irreducible representation and appears dj-timcs in the decomposition. Hence any unitary R{g) G 'R-{G) 
is decomposed as R{g) = 0^^^ Rjid) ^dj- The C*-algebra generated by TZ{G) will be denoted by Alg{TZ{G)) and 
consists of elements of the form A = J^g , G C which can be written as (12). 



C. Decoherence-free subalgebras 



For an open system with a time-independent Hamiltonian interacting with a stable environment being at the 
equilibrium state (heat bath) we expect that for an arbitrary initial state the system returns to its thermal equilibrium 
state determined by the temperature of environment. In mathematical terms it means that F* (or T*) possesses a 
unique stationary state Peq and lim„^oo r"(A) = Tr(peg^)l (liiiit_>oo TtiA) = TY{peqA)l). We say that the dynamics 
F (Tt) is ergodic in this case. The relaxation process can be very complex and involving many time scales of very 
different magnitude, some of them being much longer than the time scale relevant for the operation of our quantum 
device. Therefore, it is often very useful to study simplified approximative noncrgodic dynamics which take into 
account only the fastest relaxation processes. For such a dynamics we are interested in observables and states which 
evolve according to the reversible (unitary , Hamiltonian) evolution i.e. in a "decoherence-free" way. We shall prove 
that such observables are hermitian elements of a certain unital C*- subalgebra of Mn- In practice, the manifest 
form of this subalgebra is difficult to find and quite often we must be satisfied with a detailed knowledge of its certain 
nontrivial subalgebras. It turns out that the symmetries of the system can be very useful in this context. 

The mathematical techniques used here are based on the earlier results on the ergodic theory of quantum dynamical 
semigroups [19,20]. Some partial results of this type can be found in the recent literature also [21]. 

A minimal condition which should be satisfied by a decoherence-free observable A with respect to the dynamical 
map r is the lack of dissipation expressed by 

T{A^) = (r(A))^ , A = A* . (13) 

We formulate this condition in terms of a dissipation function Dr{., .) from Mn x Mn to Mn defined by 

DriA, B) = T{A*B) - V{A*)T{B) (14) 
Denote by A/R,(r) the following subset of hermitian matrices (observables) in A^„ and by A/'(r) its complexification 

AAR(r) = {A = A*€ Mn, Dr{A, A) = 0} , 



Af{T) = {A + iB;A,BeAfnir)} . (15) 

We shall call A/'(r) the decoherence free subalgebra of F and this name is justified by the following theorem 
Theorem 1 

1) A/'(r) is a unital C*-subalgebra of Mn 

2) F restricted to A/'(r) is a unity preserving homomorhism, i.e. for any A,B£ A/'(r) 

r{AB) = r{A)T{B) , r(i) = i (i6) 
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and moreover F is invertible on r(7V(r)). 

3) There exists a unitary Up & Mn such that 

T{A) = U^AUr , for all A e JV(r). (17) 

Proof We show first that for any A G ^R(r) and arbitrary B G A4n 

T{AB) = T{A)T{B) , T{BA) = T{B)r{A) . (18) 

The relations (18) follow from the Kadison inequality (9) applied to ^ + XB with an arbitrary A e C. Namely, 

< Dt{A + XB,A + XB) 



= r{A^) + \X\^T{B*B) + Xr{AB) + Xr{B*A) 

-T{Af - \Xfr{B*)T{B) - Xr{A)T{B) - Ar(s*)r(^) 

= \Xf{r{B*B)-T{B*)T{B)) 

+X{T{AB) - T{A)r{B)) + X{T{B*A) - T{B*)r{A)) . (19) 
Dividing both sides of the inequality (19) by |A| and taking |A| — > we have for any 2; e C, |^;| = 1 

z{T{AB) - T{A)T{B)) + z{T{B*A) - T{B*)T{A)) > (20) 

what implies (18). 

It is now a simple exercise to prove using (18) that if A,Be A/'(r), then A* ,B* , A + B, AB e A/'(r). Hence 
A/'(r) is a unital C*-subalgebra and from the general results on finite dimensional C*-algebras reviewed in Section 2.2 
it follows that (17) holds also. 

Theorem 1 gives us the general mathematical structure of observables evolving in a reversible way under a single 
step dynamics. The physical meaning of the algebra Af{T) is illustrated by the most elementary physical process 
consisting of the following stages: 

1) preparation of the initial state p, 

2) irreversible evolution governed by the dynamical map F, 

3) measurement of a given observable with a spectral resolution 

^ = flfeffc • 

k 

The probability of obtaining the value is equal to Pk = Tj:{pT{Pk)) and if ^ e A/'(r) then pk = Ti{pUpPkUr)- 
In the language of quantum information theory it means that if we choose decoherence-free decoding observables then 
the noisy quantum channel is equivalent to a certain noiseless one. 

Analogical structures can be defined for the discrete (continuous) time dynamical semigroups {F*^, A; = 0, 1,2, ...} 
{{Tt, t > 0}) to obtain the global decoherence-free subalgebras which contain observables evolving in a reversible way 
for all times 

00 

^(F(-)) = fl A/-(F'=) , ^(r(.)) = fl Af{Tt) (21) 

k=l t>0 

respectively. 

As a direct consequence of this definition, the Theorem 1 and the general results of Section 2.2 we obtain 
Collorary 7V(r(')) and A/'(T(.-)) are globaly invariant under the action of F*^ and Tt respectively. The both semi- 
groups restricted to decoherence-free algebras produce groups of unitary automorhisms {Up, k E Z} and {Ut; t G R}. 
Moreover if for the discrete case Ur can be continuosly transformed into identity then Ur S 7V(F( )). Analogical result 
holds automatically for continuous time semigroup, i.e. Ut G AA(T( '). 
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For the future applications it is useful to define a set of fixed points of T as T{T) = {A € A^„;r(j4) = A}. In 
general, ^{T) is a self-adjoint linear subspace of Mn but under additional assumptions it becomes C*-algebra. For 
the dynamical semigroup ^(T(.)) = Ker{L) = {A G Mn, L{A) = 0}. 

The algebras 7V(r), M{r'--'^) and AA(T(.)) arc the largest ones with the dccohcrcncc-frcc property but their definition 
is not constructive and often difficult to apply. Therefore, it is sometimes convenient to consider smaller subalgebras 
which are easier to construct explicitely. Theorems 2,3 provide us with the appropriate examples. 

Theorem 2 

Assume that T = UTd where U{X) = U*XU is a unitary automorphism and 

r 

Td{X) = J2w*XWc, . (22) 

Q = l 

Then the commutants defined below are C*-subalgebras and 

Wi = {W„, W*; a = 1,2, r}' C W2 = {W^W^; a, (3 = 1, 2, r}' C Af{T) . (23) 

If additionally U commutes with Tn then 

Wi C W3 = {Wc^Wp, W^W;, W*W;; a, 0=1, 2, r}' C ^(r(-)) . (24) 

For continuous time dynamical semigroups we have: 
Theorem 3 

Assume that the semigroup generator (8) decomposed as L = Lh + Ld,Lh = i[H, ■] satisfies 

LffLn = LjjLh ■ 

Then the commutant defined below is a C*-subalgebra and 

{y,-,y;;i = l,2,...,p}'c^(T(.)). 

Proof of Theorems 2 and 3 Prom the definition of the commutant Wi C W2 C M{Td) = ■A/^(r). 

Due to the condition UVd = TdU, N{T^j^) = A/'(r( )). Finally, if X = X* commutes with all W^W^ , WaW^ , 

W*W^,a,P = l,2,...,r then r''{X^) = (r'^(X))^ for aU k = 0,1,2,... 

Similarly, for continuous time semigroups the commutation of the Hamiltonian part Lh with L o implies that the 
elements of the commutant {Vj, Vjij = 1, 2, which lie in a kernel of Lu, are decohercncc-frcc. 

Remark The decomposition of the dynamical map into reversible and dissipative parts in Theorems 2,3 is not 
unique and can be optimized to obtain the largest commutants. 

To apply Theorems 2, 3 it is often convenient to use the results of the theory of group representations. Assume 
that the matrices Wa belong to Alg{TZ{G)) C A4n of a certain unitary representation of a finite or compact Lie group 
G. In this case the Hilbert space is decomposed into carriers of irreducible representations of G as in (11) where nj 
are dimensions of irreducible representations and dj their multiplicities. The commutant 

m 

{7^(G)}' = In, ® Md, C Wi (25) 

in this case. 

Similarly, for quantum dynamical semigroup satisfying the assumptions of the Theorem 3 the matrices Vj.V* ,j = 
1, 2, ...,p generate a certain Lie algebra in Mn and hence the representation TZ{G) of the corresponding group G. In 
this case the commutant (25) coincides with {Vj,VJ';j = 1, 2, ...,p}'. 

Another application of the group theory is possible when the finite or compact Lie group G is a symmetry group 
for the dynamical map T i.e. 

TZgYTZ-' = r , 7^g(X) = R{gyXR{g) (26) 

for any unitary R{g) € TZ{G), g G G. We call the condition (26) global invariance of the dynamical map T with respect 
to the group G. Using the definition of Af{T^'^) and then applying (26) we can prove that from A £ A/'(r( )) it follows 
that 'JZg{A) e A/^(r^'^) too. Therefore, the invariance holds 

7^<,(AA(^(■))) =AA(r(-)) . (27) 
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If G is a connected compact Lie group then any R{g) can be continuosly connected to identity and hence by general 
theorems from Section 2.2 TZg is an inner automorhism of ^^(r^')), i.e. R{g) G ^^(r^'^). Hence the C*-algebra 
generated by TZ{G) is included in in the decoherence-free algebra ^^{T^'^). The situation for discrete symmetry may 
be quite different due to "incidental degeneracies" of F which admit larger group algebras. All that is true for the 
dynamical semigroup {Tf, < > 0} also and will be illustrated by examples in Section 3. 

The stronger invariance conditions with respect to the symmetry group called local invariance can be formulated 
for the dynamical map F = UTn = as in the Theorem 2. It reads for all a and g G G 

R{g)*Wo,R{g) = (28) 

Therefore the algebra generated by the group representation Alg{TZ{G)) C J^{r'^'^). Similarly, for quantum dynamical 
semigroup as in the Theorem 3 the condition 

R{g)*Vo,R{g) = (29) 

implies Alg{H{G)) C A^(T(.)) = KeriLo). 

Summarizing, global invariance (26) for connected compact Lie group implies Alg{n{G)) C 7V(r( )) while local 
invariance (28) or (29) implies for any compact group G that Alg{TZ{G)) C 7V(F(')) or Alg{7l{G)) C ^f{T(^.)) respec- 
tively. 

The concrete examples of semigroups and their decoherence-free observables obtained using the methods proposed 
above will be discussed in Section 3. 

D. Limited relcLxation 

We expect that the open system governed by a non-ergodic dynamics shows a behaviour which can be called limited 
relaxation. It means that in the Heisenberg picture any initial observable becomes for long times decoherence-free and 
the set of decoherence-free observables is nontrivial. We are able to prove this property under additional assumptions 
which are typically fulfilled for open quantum systems weakly coupled to heat baths. 

We begin with the dynamical map F satisfying the following conditions: 

1) There exist a faithful (i.e. given by a strictly positive density matrix) stationary state a , F*(o") = cr. 

2) Detailed balance condition with respect to a holds, i.e. F = UTd = TdU where U{X) = U*XU is a unitary 
automorphism ,h({a) = a and the dynamical map F/j is a hermitian operator on the Hilbert space {Mn \ < , • >cr) 
equipped with the scalar product 

<A,B >„= Ti{(jA*B) , A,B€Mn (30) 

and the norm |m|<T =< A, A >^/^ . This Hilbert space is often called Liouville space. More explicitely U is unitary, 

T D is hermitian and IATd = T]jU is a normal contraction on the Liouville space(A^„; < •, • >ct)- 

The existence of a stationary state , generally non unique, is quaranteed by the finite-dimensionality of the system's 

Hilbert space but its faithfulness is a nontrivial assumption. 
Lemma 1 

If the dynamical map F possesses a faithful stationary state a then the fixed point set JF(F) is a C*-subalgebra of 
Mn and jr(F) c^(F( )). 

Proof Take A = A*e ^(F) then, A"^ = (F(A))^ and hence 

= Tr{a[A^ - F(^)F(^)]) = Tr{a{T{A^) - F(A)F(A))) . (31) 

As r{A^) - T{A)T{A) > and a is faithful then (31) is equivalent to r(A2) = {TiA)f . It follows that A e Af{T). If 
A = A* € A/'(F) then we can reverse the reasoning to prove that A € ^(F). For a non- hermitian A we can use the 
decomposition A = Ai + iA2 into hermitian ones. 

We shall denote by Pr the orthogonal projector acting on the Liouville space {Mn, < •> • >cr) and projecting on the 
subspace Af{r^'^). 

Limited relaxation process is characterized by the following theorem. 

Theorem 4 

Consider a dynamical map F satisfying the detailed balance condition of above. Then M{T^'^) = J^{T]j) and 

lim ||F*^(A) - U''{PrA)\\ = , for all A G Mn ■ (32) 

k—>oo 
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Proof From T'' = U^V\, it follows that AA(r( )) = N{rj^). For any A G 7V(r( )) we have ||r^(A)||<, = 

1 1 j4 1 1 Therefore, because To is a hermitian contraction on the Liouvillc space, M{T^j^) = J^{Td)- 

For the dynamical semigroup with the generator L = i[H,-] + Ld the detailed balance condition means that 
the Hamiltonian part Lh commutes with L]j, [H,<t] = and is a hermitian operator on the Liouvillc space 
{A4n', < ■, ■ >a) [7]. Introducing a kernel Ker(L£)) = {A £ A^„;L£)(A) = 0}, the orthogonal projection on Ker(L£)) 
denoted by Plo *ind applying the same idea of the proof we obtain the following analog of the Theorem 4. 
Theorem 5 

Consider the dynamical semigroup Tt with the generator satisfying detailed balance condition of above. Then 
A/'(r(.)) = Ker(L£,) and 

lim ||rt(^) - e''"PL„{A)e-'*"\\ = , for all A G Mn ■ 

t — ^oo 

Example 

The particular case of detailed balance generator with the stationary Gibbs state 

a = Z-^e-"/^ (33) 

leads to a very specific structure oi Ld [7] 

Lo{A) = (^/^^^ - \{y;y3^A} + e-^/^{V,AV* - ^-{V,V;,A})) (34) 

with [H, Vj] = u)jVj . Simple calculation involving the dissipativity condition (10) and hermicity of Ld on the Liouvillc 
space shows that in this case 

< Tr(CT[L(A*A) - L{A*)A - A*L{A)]) = - < A, Ld{A) 

= ^ (< [y,-. A], [V^,A] +e-^/^ < [V*,A], [V*, A] . (35) 

j;wj>0 

Therefore, the commutant {Vj,V*;ujj > 0}' = Ker(LD) = A/'(T(.)) . 

III. EXAMPLES 

In this Section we give examples of physical systems governed by nonergodic quantum dynamical semigroups and 
illustrate the different methods of searching for decoherence-free observables introduced in Section 2. 

A. N-particle systems with permutation invar iance 

Consider a system which consists of N identical particles each of them is described by a finite dimensional Hilbert 
space H^^^ = C^; j = 1, 2, N. There are two possibilities of constructing a quantum dynamical semigroup generator 
invariant with respect to the permutation group Sn- The first one corresponds to particles interacting with identical 
"private" reservoirs and is given by the Master equation in the Schrodinger picture for the A/'-particle density matrix 
Pt 

d ^ 

-p^ = -i[H,pt]+Y,L*M (36) 

m=l 

where H isaN- particle permutation invariant Hamiltonian and LJ^ are identical copies of a single-particle dynamical 

semigroup generator. 

The second example corresponds to the collective coupling of iV-particle system to a single "common" reservoir 

j^pt = -i[H,pt] + \ ^([Fa,ptV„] + [Vo^pM) (37) 

a 
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where Va = X]rn=i ''a with identical copies of a single particle operator Va- 

The Hilbert space of the system possesses the structure of A^-fold tensor product with the basis {eji,j2,---jN = 

Sji <8) (8) • • • (8) Cjjv}. Any permutation tt G Sn is represented by the operator -R(7r) defined in terms of the basis as 

■R(^)ejlj2,...J.Y = e7rOi),7r02),...7rO_v) • (38) 

It is well-known that this representation is highly reducible and therefore Alg{TZ{SN)) is nontrivial. The explicite 
construction of the decomposition (11) (12) is given in terms of Young tables and described in the texbooks on group 
theory [22]. 

The semigroup generated by the master equation(36) is globally invariant with respect to S'jv but because S'jv is 
a discrete group according to the results of Section 2.3 we can have no nontrivial decoherence-free observables i.e. 
A/'(T(.)) = C. Indeed, if a single-particle semigroup generated by is ergodic and the Hamiltonian H = 0, then the 
whole A^-particle dynamics governed by (36) is ergodic too. 

For the case of the dynamics goverened by the master equation (37) the situation is different. The corresponding 
semigroup is locally invariant with respect to Sn and therefore according to the results of Section 2.3 Alg{TZ{SN)) C 

A/-(r(.)). 

B. Superradiance model with SU{2) symmetry 

The superradiance [23] and subradiancc phenomena [24] are the first examples of collective dissipative phenomena in 
quantum systems which leads to a nonergodic behaviour (on a certain time scale) due to (approximative) permutation 
symmetry of the interaction of atoms with the electromagnetic field. In the simplest case of 2-level atoms in open 
space and for the temperature T = 0, we obtain the master equation of the form (37) with 

1 ^ 1 
H=-ojY1 ^i™^ ' ^« = 2^i^^ + '^v) (39) 

m=l 

where w is an atomic frequency, 7 is a radiation damping constant and <j'j^\j = x,y,z, are Pauli matrices for the 
m-th atom. 

In order to find decoherence-free observables we can use cither the local symmetry of the master equation with 
respect to the permutation group Sn or the fact that the collective operators J2m=i ''"j "'' which appear in the 
semigroup generator define a reducible A''-fold product representation of the group SU{2) on 0^ C^.Then we obtain 

AlgiTZiSN)) C Alg{n{SUi2))) C AA(T(.)) (40) 

where TZ{Sn) is a standard representation of Sn defined by (38) and TZ{SU{2)) is the mentiond above product 
representation of SU{2). For our case of T = the exphcite structure of A/'(T(.)) is unknown while for T > we can 
use detailed balance condition and the Theorem 5 to characterize A/'(T(.)) [20]. 

IV. NON-MARKOVIAN CONTROLLED OPEN SYSTEMS 

The analysis of the previous Section has a rather mathematical and phenomenological character and its most 
developed part is applicable to dynamical semigroups (discrete or continuous) and hence involve Markovian assumption 
or at least stationary external conditions. For a generic controlled system we expect that the external conditions like, 
for instance, external electromagnetic fields used to control the quantum device vary in time. Therefore, we need other 
approximation schemes which are based on the assumption of the weak influence of the environment on the system 
and express the dynamics in terms of the fundamental ingredients (controlled Hamiltonian, correlation functions of 
the bath) [25]. 

A. Errors in CQOP 

We would like to describe a quantum device operating during the time interval [0, t]. We assume that in the initial 
moment the prepared state can be given by a product formula (1) with a pure state p = JV" >< while in the final 
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moment the state of our open system is described by the reduced density matrix (2) . The actual dynamics is compared 
with the ideal one given in terms of the controlled unitary evolution Us{s) which solves the Schrodinger equation 
governed by the controlled Hamiltonian Hs{s); < s < t. The error e due to the interaction with environment can 
be defined as 

e = l-<Usm,At{\i^><i>\)Usm> . (41) 

The choice of the Hamiltonian Hs (t) is a tricky problem. Typically, one takes a sum of the time-independent Hamil- 
tonian Hg of the bare system S and a time-dependent contributions Hg{s) describing the influence of controlled 
external fields on our system. However the interaction with environment induces generally time-dependent Hamil- 
tonian corrections. For example, vacuum fluctuations of the electromagnetic field produce Lamb shift of the atomic 
energy levels and van der Waals type interactions between neutral atoms [24]. Similar phenomena appear in solid 
state physics due to interactions with phonons and other quasi-particles. In the following we adopt the optimistic 
point of view assuming that we can include these corrections in the process of designing quantum algorithms. In 
practice one should expect the presence of the Hamiltonian errors due to an approximative knowledge of the true 
physical Hamiltonian Hs{s). 

B. Reduced dynamics in Born approximation 

We consider the system which consists of the quantum device S driven by the time-dependent Hamiltonian Hs{t) 
and the bath with its Hamiltonian Hn. Then the total Hamiltonian reads 

H{t)=Hs{t)+HR + Hi^t (42) 

with the interaction Hamiltonian of the form 

Hint =^So.(^Rcc (43) 
a 

where Sa,Ra arc self-adjoint operators. 

In the following we investigate the dynamical map F* in the Schrodinger picture which describes the evolution of 
S from the initial moment tin = —t to to the flnal one tfin = r in Born approximation. The exact dynamical map 
can be expressed in terms of Hint, the unitary propagator for the total systems U{t,s), the free propagator for the 
total system Uo{t, s) and the propagator Us{t, s) defined by 

U{t, s) = Te-'Il ^'"^ , Uo{t, s) = Te-'/^^^("'+^«) '^^ , 

C/s(i,s) = Te-'./'>^(")''" (44) 

where T is the time ordering operator. Introducing the "supcropcrator' notation: Up= UpU* and HintP = [Hint, p], 
choosing the initial state similarly to (1) 

Psr{-t) = p(E)ojR, [Hr, ojr] = 0, Tr{ojRRa) = (45) 

and applying the second order integral identity for U (i, s) we obtain the exact equation 

r*{p)=Us{T,-T)p 



-Tr|y^ ds j duUo{T,s)HintUo{s,u)HintU{u,-T)p<S'COR^ . 



(46) 



The Born approximation consists in replacing U{u, —t)p®ujr by Us{u, —t)p®lor in eq.(46). Our main assumption is 
that for a given initial state p and a given trajectory {Hs{t); — r < t < t} of the controlling Hamiltonian the influence 
of the reservoir on the state of our system S given by the second term on the RHS of eq.(46) is small, say of the order 
e << 1. Then up to the higher order terms ^ we can use the proposed approximation. One should notice that we 
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do not necessarily impose weak coupling or/and short time r regime but we can search for the optimal choice of p 
and Hs{t) to fulfill the Born approximation beyond this regime. 
The final formula for T* (p) in Born approximation reads 

V*{p) = Us{p + ^*{p)-\{K,p}) (47) 

where tJs = Us{t,—t). The superoperator called error map is a completely positive map given by 

^*{p) = J dsj duTv{ujRRaR0{s-u))Sf3{s,-T)pSc,{u,-T) (48) 

where Sa{s,u) = Usiu,s)Sa , Ra{t) = e^^^'^Ra- The operator isT > is given by if = $(1). Strictly speaking one 
has also a contribution of the form — p] which contains Hamiltonian corrections due to the interaction with R like 
for example Lamb shift, collective Lamb shift, etc. but we put h = Q. This can be justified by starting with a bare 
Hamiltonian Hg{t) containing appropriate counterterms and performing renormalization procedure. It is , however, 
equivalent to treating Hs{t) as a full physical Hamiltonian and puting h = 0. 

The equation (47) resembles in its structure the short time approximation to quantum dynamical semigroup and 
indeed, if the conditions of Markovian approximation are fulfilled we can derive from (47) different types of the 
Markovian master equations [7,9,11]. The approximative map (47) is trace preserving but complete positivity is 
satisfied up to the higher order terms in e. 

Finally, we pass to the frequency (energy) domain introducing spectral density of the reservoir Rai}{ijj) by 

/oo 
Raf3{u^)e-^'(L; (49) 
-oo 

and defining 

Y^{uj) = J So.{s,-T)e-''^'ds . (50) 
We obtain the following formula for the error map 

/oo 
du;R„0{uj)Yf,{uj)pY:{u;) (51) 

which is a kind of quantum fluctuation-dissipation theorem in the linear response regime [9,10]. 

C. Error formula 

The general definition of the error due to the interaction with the environment (41) can be applied to our approxi- 
mative evolution 

e = 1- < Usip, r*(|V >< ipDUsip > . (52) 
Putting the expression (47) we obtain the error formula in Born approximation 

e=<'ip,KiP> - <iP,<i>*i\^j><'iP\)'iP> (53) 

which can be the starting point for the analysis of efficiency of quantum devices. Writing it in terms of (49-51) we 
obtain 



/oo 
dtoR^^ico) (<^I;,Y:{uj)Y^{uj)^> 

- < V', Y*{u;)^ X V, F/3(a;)V >) • (54) 
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The interpretation of (54) becomes simple and elegant if we introduce the following notation 

< V.l^aMF/sMV > - < X i^,Yi3{uj)i, >= 2r5„^(a;) . (55) 

Now we can write the error formula in a compact form as an overlap of two "correlators" 'lZ{u)) = [Ra0{co)] and 
<S(a;) = [5^/3(0;)] multiplied by the duration of the operation time of the quantum device 2t 

/oo 
cL;Tr(7^(w)5(w)) . (56) 
-00 

Both correlators encode the information about the dynamics and the initial state of the baths and the quantum 
device, respectively. Moreover, one can see that their mathematical structure is essentially the same. To show this 
we consider a sequence of random (quantum or classical) time-dependent real variables fa{t) and denote by < • >f 
the appropriate average. Define a correlator J^{ui) = [Fa0{w)] by the following formula 

^lim_ — J ds(< Uit + s)fpis) >F-< fa{t + s) >f< fp^s) >f) 

/OO 
F„/3(w)e-'"*(iw . (57) 
-00 

Then puting /„ = < • >f= Tr(a;/{-) we obtain = 'R,{uj). Similarly, T{uj) = S{uj) for /„ = Sa, < ■ >f=< 

il), -tj-j > if only T is long enough such that the limit t ^ 00 in (57) can be used. 

Having the formula (56) one can try to design different strategies reducing errors. Although very roughly the error 
seems to be proportional to the operation time 2r the situation is much more complicated. For example, if the leading 
error source is described by the correlator satisfying for small uj 

n{u)) ~ w'" , K>1 (58) 

the optimal strategy is to slow down the time evolution of S in order to move the correlator S{u)) into the low frequency 

domain. On the contrary, if TZ{uj) is flat "fast gates" are optimal. Generally, one can try to shift the support of S{u!) 
into the frequency domain where TZ{uj) is small. However, because both correlators are positively defined matrices 
we cannot expect dramatic reduction of errors by some "cancelation effects". Moreover, different type of reservoirs 
and interactions produce different shapes of 'R,{lj), e.g. scattering with gas particles gives flat shapes while the linear 
coupling to quantum bosonic field produces the relations of type (58) with some cut-offs for large uj (see [25] and the 
authors contribution to [11]). 

Another different strategy is to select an optimal subspace of initial states ip which are " decoherence-free" with 
respect to the leading error source. The following simple criterion is a generalization of the idea of quantum correcting 
codes [12]. 

Assume that the matrics 'lZ{oj) is strictly positive for a; e f2. Than any ij) satisfying the eigenvector condition for 

all a and w G 

F«(w)V' = Aa(w)V' , Aa(w) e C (59) 

evolves without errors. 

Another partial results and examples of applications can be found in [25] . 

V. CONCLUSIONS 

The author believes that the presented results, although often formulated in a rather abstract mathematical form, 
can serve as a guide in the discussion concerning the optimal strategies in controlling quantum devices. These strategies 
can be roughly divided into passive and active ones. At present, it seems that the passive ones can be realized in 
three apparently different ways: 

I) using the symmetries of the coupling to the baths and systems' Hamiltonian to obtain large enough decoherence- 
free subalgebras as presented in Section 2.3, 

II) slowing down the action of quantum algorithm (slow gates) for the environments which satisfy the low frequency 
scaling (58), 
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Ill) Speeding up the quantum algorithm (fast gates) for the environments with flat 71{lu) ~ 7?.. 

The best known models of the type I) are given by the " superradiance" systems governed by the equations (37). In 
principle, they can be realized in different settings using collective coupling to bosonic fields (photons, phonons, etc.,). 
Such reservoirs satisfy also the condition (58) what suggests that both I) and II) demand the similar technology. 
Indeed, the collective coupling of qubits to the bosonic field, which is a necessary condition for the validity of (37) is 
possible if the wave-length of the relevant bosonic modes is long in comparison with the diameter of the device [24]. 
It means that only the low frequences can appear in the correlator S{uj) (55), (56). This is equivalent to the slow 
gates strategy II). Moreover, in both cases I) and II) the time scale of gates must grow with the size of the device 
(proportional to the number of qubits) in order to keep the total error constant. 

The main obstacle for the strategy I) is the approximative character of the symmetry which becomes less accurate 
with the growing size of the device. Moreover, the Hamiltonian corrections due to the renormalization effects discussed 
in Section 4.2 do not possess this symmetry except of the cases with the very special geometry [24]. 

The strategics II) and III) arc obviously contradictory and therefore for real systems , for which both types of 
reservoirs are at work, the optimization is needed. 

The most elaborated active strategy is the idea of quantum error corrections for quantum computers [12]. The 
existing models of error corrections possess, however, several drawbacks. First of all, the noise is assumed to be 
independent of the algorithm what is equivalent to the assumption of the strategy III) (7?.(w) — TZ)- Moreover, the 
entropy produced in the computer must be very precisely taken away by "fresh" qubits prepared in pure states. 
Finally, initial errors for all qubits used in this scheme are not taken into account seriously. 

Another active approach, called dynamical decoupling [13], involves the existence of the frequency domain for which 
the correlator TZ{lu) ~ 0. Then applying a suitable fast modulation of the system's dynamics we can, in principle, 
move the support of the correlator S{lli) into this domain. Although for the specific environments the correlator TZ{ijj) 
may have some "valleys" , the overall tendency for Ti{uj) is to grow with u. This is related to the fact that the number 
of accesible quantum states of the environment typically grows with increasing energy scale. In any case, dynamical 
decoupling strategy demands fast gates technology as for III). 

In authors' opinion the optimized combination of the passive strategies I, II and III may be the most realistic 
attempt to deal with quantum noise and associated errors. It would be interesting to apply the formulas of the 
type (56) to estimate the minimal error per a single gate within the different proposed physical implementations of 
quantum computers (e.g. trapped ions , NMR , quantum dots , etc). 
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